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ABSTRACT 

Recent observations in galaxies and clusters indicate dark matter density 
profiles exhibit core-like structures which contradict to the numerical simulation 
results of collisionless cold dark matter. On the other hand, it has been shown 
that cold dark matter particles interacting through a Yukawa potential could 
naturally explain the cores in dwarf galaxies. In this article, I use the Yukawa 
Potential interacting dark matter model to derive two simple scaling relations on 
the galactic and cluster scales respectively, which give excellent agreements with 
observations. Also, in our model, the masses of the force carrier and dark matter 
particle can be constrained by the observational data. 

Subject headings: dark matter, galaxies, clusters 
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1. Introduction 

The nature of dark matter remains a fundamental problem in astrophysics and cos- 
mology. The rotation curves of galaxies and the derived mass profiles in clusters indicate 
the existence of dark matter. It is commonly believed that dark matter is collisionless and 
becomes non-relativistic after decoupling. Therefore, they are regarded as cold dark mat- 
ter (CDM). The CDM model can provide excellent fits on large scale structure observations 
such as Lyg spectrum ( Croft et al.lll999l : ISpergel and Steinhardtll2000[). 2dF Galaxy Re dshift 
Survey ( jPeacock et al.l 120011 ) and Cosmic Microwave Background (ISpergel et al.l 120071 ). 



However, N-body simulations based on the CDM theory predict that the density profile 
of the collisionless dark matter halo should be singular at the center (p y ?"~^ ') (iNavarro et al. 



19971 ) while observations in dwarf galaxies give p 



n~j f 



-0.29±0.07 



(lOh et al.ll201ll : lLoeb and Weiner 
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201ll ). On the cluster sc ale, observational data from gravitationa l lensing also sho\ y that cores 



exist in some clusters (JTyson et al.l Il998l : iNewman et al.ll201ll ). In particular, ISand et al. 
( 120081 ) get a = —0.45 ± 0.2 by the combination of gravitational lensing and dynamical data 
of clusters MS2137-23 and Abell 383. Clearly, observations do not support the numerical 
small-sca le predictions by the CDM model. This discrepancy is known as the core-cusp 
problem Jde Bloklboipf ). 



Many theories have been invoked to solve the core-cusp problem. For example, some 
bary onic processes such as supernova feedba ck have been suggested to alleviate the prob- 
lem ( iWeinberg and Katzll2002l : Ide Blokll2010l ). Recent high resolution cosmological hydro- 



dynamical simulations show that cored dark matter density profile could be produced in 
Milky- Way like halo i f there is enough ra diation pressure of massive stars before they ex- 
plode as supernovae (JMaccio et al.l 1201 ll ). However, it is still controversial to make the 
conclusion because t he actual contribution o f supernova explosions is limited by the low star 
formation efficiency (jPenarrubia et al.l 120121 ). Also, it is challenging to i nvoke baryonic pro 



cesse s as the main mechanisms for some dark matter dominated galaxies (jVogelsberger et al. 



20121 ). Another spectacular idea is that dark matter is not cold. The existence of keV sterile 
neutrinos, as a candidate of warm dark matter (WDM)^ has been proposed to solve the 
discrepancies (JBode et al.ll2001t IXue and Wij|200ll : IChol l2012l ). However, recent observa- 
tions tend to reject the keV sterile neutrinos to be the major component of dark matter 
since the observational bou nd of sterile neutrino mass in Lyraan-alpha forest contradicts to 



that in x-ray background (JAbazajian and Koushiappad l2006l : IViel et al.l l2006l : ISeljak et al. 



20061 ) . Als o, the WDM model alone cannot get a good agreera ent on the large scale power 
spectrum (ISpergel and Steinhardtl l2000l : iBoyarsky et al.l 120091 ) . Therefore, the success of 
the CDM model on large scales suggests that a modificati on of the dark matter proper - 
ties may be the only approa ch to solve the discrepancies (ISpergel and Steinhardtl I2OOOI ) . 
Spergel and Steinhardtl (I2OOOI ) proposed that the confiict of ob servations and s imulations can 



be rec onciled if the CDM particles are self- interacting. Later, iBurkertI (J2000l ): lYoshida et al. 



(I2OOOI ) performed numerical simulations of self-interacting dark matter (SIDM) with con- 
stant cross section and showed that core-like structures could be produced. However, this 
proposal fell out of favour because gravitational lensing and X-ray data indicate that the 
cores of cl usters are dense and elli psoidal where SIDM model predicts that to be shallow and 
spherical ( jLoeb and Weinerll201ll ). Also, there is a discrepancy on the required cross section 
per u nit mass in dwarf galaxies (a/m ~ 04.~ -^ '^P^^ ^~^^ ^^^ clusters (a/rn ~ 0.01 — 0.1 cm^ 



g-^) flMiralda-Escudel I2OO2I : iRandall et al.ll2008l : iBucklev and Foxll2009l : iTulin et al.l 120131 ). 
Recent numerical simulations indicate that only a small wind ow open for a consta nt cross sec- 



tion SIDM model to work as a distinct alternative to CDM (jZavala et al.ll2013l ). Therefore, 
a velocity-dependent cross section of the SIDM were explored to tackle the problems. 
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Loeb and Weinerl (120111 ) proposed that the possible existence of a Yukawa potential 
among the dark matter particles can resolve the problem raised by the SIDM with constant 
cross ssection. The velocity dependence could naake s cattering important in dwarf galaxies 
but unimportant in clusters (jVogelsberger et al.ll2012[ ). In this model, the dark matter par- 
ticles of mass m is set by an attractive Yukawa po tential with couphng strength a mediate d 



by a gauge boson of mass ra^ in the dark sector ( JFeng et al.ll201 
The cross section is well fitted by 

^/32ln(l + r'), /3<0.1, 



Loeb and Weineiil201ll ). 



a 



Svr 



■/3V(1 + l-5/3^-^^), 0.1</3<103, 



-^(ln/3 + l-0.51n-^/3)2, /3 > 10^, 

where (i = 7iv^„^/v'^ = 2 am^/{mv'^) and v is the relative velocity of the dark matter particles 
(JLoeb and Weinerl 1201 ll ). The Vmax is the velocity at which av reaches its maximum value. 
In the original proposal of the Yukawa interacting dark matter model, the parameter Vmax 
was considered to be v^ax ~ 10 — 100 km/s in order to match the characteristic velocities of 
dwarf galaixes and clusters so that the scattering rate is larger in dwarf galaxies but much 
smaller in clusters. Nevertheless, this is not a necessary range. In Fig. 1, we can see that 
the cross section per unit mass is still larger on the galactic scale (f ~ 10 — 100 km/s) even 
if Vmax ~ 10^ km/s. Also, the scattering rate of dark matter particle is pcio lrn)v, which is 
still larger in galaxies because the central density p^ is about 10 times larger than that in 
clusters. In this article, I release the free parameter v^ax to ~ 10^ km/s so that the cross 
section will be cr oc (ln/3 + 1)^ and a oc /3~°-^^ on the galactic (/3 > 10^) and cluster scales 
(/3 ~ 10^) respectively. Therefore, the cross section is more velocity dependent on the cluster 
scale. In the following, I will use the above velocity-dependent self-interacting dark matter 
(vdSIDM) model with Vmax ~ 10^ km/s to derive a scaling relation to relate the total mass of 
dark matter halo M and v. Then, I will compare the derived scaling relations with empirical 
fits on the galactic and cluster scales. 



2. vdSIDM model in galaxies and clusters 



In the vdSIDM model, the size of a core Re in a structure depends on the self-interacting 
rate of dark matter particles. Inside the core, we may assume that the dark matter particles 
interact with each other at least once during the evolution of a galaxy or cluster. Therefore, 
at r = i?c, we have 

p, (-) V ^ i/o, (2) 
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where pc = 3Mc/47ri?|? is th e central density of the core, Mc is the total mass of the core and 
Hq is the Hubble constant (JRocha et al.l 120121 ) . By using the Virial relation v 
and Eq. (2), we get 



'"'-i^y\B"""^-"'^-"-'^"- 



VGMJR, 
(3) 



In ga laxies, since the core mass is about one-tenth of the total mass {Mc 
2012[ ). we have 



O.IM) (IRocha et al 



log 



M 



3.65 + 0.5 log 



a/m 



+ 3.5 log 



(4) 



For P > 10^, the cross section in Eq. 
V (see Fig. 2): 

log (-)= -0.37 log 



0.1 cm^g^iy ' ^\^lkm/s^ 

1) can also be approximated by using a power law of 



+ log(^ 
m 



(5) 



1 km s^^ 

where ctq is a constant which depends on m,^. By putting the above equation into Eq. (4) 

we get 

M \ / o'o/'m 



log 



M, 



Q 



3.65 + 0.5 log 



0.1 cm^g 



+ 3.3 log 



1 km/s 



(6) 



On the other hand, the empirical fit of the baryonic Faber- Jackson relation obtained from a 
represe ntative sample of 436 galaxies is given by log(v/l km s~^) = 0.299 log(MB/M0) — 
1.053 (jCatinella et al.l 120121 ). Since Mb ~ 0.17M, the observed Faber- Jackson relation 
becomes log(M/M0) = 4.29 + 3.341og(f /I km s^^). Compare the empirical fit with Eq. (6), 
we get (Jo/m ^ 1.9 cm^ g^^. By Eq. (5), the cross section per unit mass for dwarf galaxies 
(f ~ 50 km/s) and Milky-Way size galaxies {v ~ 200 km/s) are a/m ~ 0.4 cm^ g~"^ and 
a/m ~ 0.2 cm^ g~^ respectively. Therefore, both the power-law dependence of the scaling 
relation {^ 3.3) and the order of magnitude of the cross s ection per unit mass (a/m ^ 0.1 — 1 



cm 



2012 



g ~^) are generally agree with the recent observations (JBuckley and Foxll2009l : iPeter et al. 
Tulin et al.ll2013l ). Futhermore, since /3 > 10^ in galaxies, by using ao/m, f» 1.9 cm^ g~^ 
and Eq. (1), we have m,'^m, ~ 7r(ln/3 + l)^/(cr/m) ~ 0.3 GeV^. By combinin g the above result 



with the estimated lower bound derived from dwarf galaxies m^ > 40 MeV (JBuckley and Fox 
2009h . we have m < 200 GeV. 



Similarly, we can apply the same model to clusters. However, since v ~ 10^ km/s 
(1 ^ /3 < 10^) in clusters, by Eq. (1), the cross section drops faster with velocity (see 
Fig. 1): 

a _ 8n /32 

m m?/m (1 + 1.5/ 

-1 



1.65\ 



m 



Vn 



V 



0.7 



(7) 



By using the result m,'i,m ~ 0.3 GeV'^, we have ctq/tti = 0.012 cm^ g ^. Therefore, the 
cross section per unit mass in clusters is a/m, ~ 0.06 cm^ g~^, which is consistent with 



-5- 



the recent observed bounds a/m < 0.1 cra^ g ^ (IMira 



Buckley and Foxl l2009l : IPeter et all 120121 : iTulin et al 



da-Escudell2002l : iRandall et al.ll2008l : 



20131 ). Since t he size of a cluster is 



about 100 times o f the core size, which is equivalent to M ~ lOOM^ ( lArabadjis et al.l 12002 
Rocha et al.ll2012[ ). by putting Eq. (7) into Eq. (3), we have 



M = im[£j\,'''G--''-v^^^[^ 



1/2 



,3.15 



Since v ~ ^/SkTjm'g^ where T and vfig are the temperature and mean mass of a hot gas 



particle, we can obtain a scaling relation 



M f^ 1.7 X lO^^M 







T 



2keV 



1.58 



(9) 



Surprisingly, this derived scaling relation gives excellent agreements with both the power 
dependence and proportionality co nstant of the empirical fits from 118 clusters M = (1.56 ± 
0.01) X 10^*Mq{T/2 keV)i-^^±°-°6 dVentimidia et al.ll2012[ l. 



3. Discussion 



The original purpose of suggesting the vdSIDM model i s to explain the observe d cores 
in dwarf galaxies without affecting the dynamics in clusters f lLoeb and Weinerll201ll ). They 
assume VmnT ~ 10 — 100 km/ s so that the maximum cross section lies on the galactic 



scale ( IVogelsberger et al.l 120121 ). In fact, this is a free pa rameter whic h depe nds on m<^ 
and m, and it is not necessary to be about 10-100 km/s. iRocha et al.l ( 120121 ) show that 
Pc ~ O.OISMq pc~^(t'/100 km s~^)~''-^^ by simulations, which means the central density of 
dark matter is higher in dwarf galaxies. As a result, the scattering rate of dark matter particle 
~ Pc{cr/m)v is always larger in dwarf galaxies than clusters even if Vmax ~ 10^ km/s. On the 
other hand, the circular velocity of a dwarf galaxy can be obtained by substituting Eq. (5) 
into Eq. (2), which gives v ^ {'inGHo/3)^-^^R^-^^{ao/m)-^-^^ . For R = 250 pc and R = 500 
pc, we get f = 13 km/s and f = 22 km/s respectively. Th ese values are consistent with 



the o bserved circular velocities on the dwarf spheroidal scale (jWalker et al.ll2009l : IWolf et al. 
20101). Therefore, th e large Vmnx can still solve the t oo big to fail problem suggested by 
Bovlan-Kolchin et al.l J201l[ ): IVogelsbereer et allJ2012[ ). 



In this article, I show that if f^ 



lO'' km/s, the cross section goes like 



,-0.37 



and ~ f ~°'^ in galactic and cluster scales respectively. The derived scaling relation on th e 
galactic scale is M ex v^'^, which agrees with observations M ex v^'^^ (jCatinella et al.ll2012l ). 
The cross section per unit mass contrained by this model is a/m ~ 0.2 — 0.4 cm^ g~^ for v = 



-6- 



50—200 km/s, which is also consistent with the observed bou nds in dwarf and normal galaxies 
(JBuckley and Foxll2009l : iPeter et al.ll2012l : iTulin et al.ll2013[ ). By applying the same model in 
clusters, we get M oc T^'^*, which again gives excellent agreements with ob servations in both 
prop ortionality constant (~ IO^^Mq) and power dependence (1.57 ±0.06) ( iVentimiglia et al. 
2OI2I ). These results provide evidences on the non-power-law velocity dependent cross section 
of self-interacting dark matter, livmax is 10-100 km/s, the derived scaling relations would be 
M oc f '^•^^ and M oc T°-^^ on the galactic and cluster scales respectively. Obviously, they do 
not match the empirical fits from observational data. Futhermore, in my model, it predicts 
m\m ~ 0.3 GeV'^. If m^ > 40 MeV, then m < 200 GeV, which is a testable range in the 
future large hadron collision experiments. 



Recently, IVogelsberger and Zavalal (120131 ) study the impact of self-interacting dark mat- 
ter on the velocity distribution of dark matter haloes and the anticipated direct detection 
signals. They find that all SIDM and vdSIDM models show departure from the velocity 
distribution of the CDM model in the center of the Milky Way halo. Therefore, different 
SIDM scenarios, including my model, might b e distinguished from each othe r through the 
details of direct detection signals in the future (IVogelsberger and Zavalal 120131 ). 
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